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1. Introduction. Subexponential distributions (Chistyakov, 1964), 
as the name suggests, have tails which decay at a subexponential 
rate, and hence, provide good models for heavy-tailed data. These 
distributions are defined by a property on their rr-fold convolutions 
which make them useful in various problems involving sums of 
independent random variables. To be more precise, we need to 
introduce some notation. For a distribution F we write F = 1 — F 
its tail function and E*” its n-fold convolution. Next if / and g are 
two functions, we write f ^ g to signify that the ratio f /g tends to 
1 at infinity. A distribution function F is subexponential if and only 
if ~ 2F, or, equivalently (Chistyakov, 1964), if F*^ ~ nF for 
all positive integers n. Thus, these distributions have particularly 
nice analytic features with respect to convolution, and, as such, 
form an excellent class that allows both generality and tractibility 
in asymptotic analyses of several stochastic models. 

When dealing with problems involving extreme values, subexpo¬ 
nential distributions have two possible limiting extremal behavior; 
for instance the Pareto type belong to the domain of max-attraction 
of the Frechet distribution, while the Weibull type belong to the do¬ 
main of max-attraction of the Gumbel distribution. Hence, as far as 
extremal behavior is concerned, some subexponential distributions 
are classified as light tail, in the sense that their extremes have a 
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limiting behavior which is of the same type as that for the extremes 
of a normal distribution say. 

In this paper, we derive tail area expansions for infinite order 
weighted sums of light tailed, yet subexponential, random variables. 
For Pareto type distribution, such expansions were obtained in Barbe 
and McCormick (2005) but their form and the proofs differ sharply 
from those in the current paper. 

Examples of applications of subexponentiality include transient 
renewal theory (Teugels, 1975; Embrechts and Goldie, 1982), ran¬ 
dom walks (Griibel, 1985; Veraverbeke, 1977), branching processes 
(Ghistyakov, 1964; Athreya and Ney, 1972; Chover, Ney, and 
Waigner, 1972), queueing theory (Pakes, 1975), shot noise (Lebe¬ 
dev, 2002), infinite divisibility (Embrechts, Goldie, and Veraver¬ 
beke, 1979), ruin theory (Asmussen, 1997; Goldie and Kliippelberg, 
1998; Tang and Tsitsiashvili, 2003), compound sums (Gline, 1987; 
Embrechts, 1985; Griibel, 1987), insurance risk theory (Embrechts, 
Kliippelberg and Mikosch, 1997), and heavy-tailed linear processes 
(Rootzen, 1986; Davis and Resnick, 1988; Geluk and De Vries, 2004; 
Ghen, Ng, Tang, 2005). In these papers, inputs to the processes fol¬ 
low a subexponential distribution and an asymptotic analysis of an 
output is obtained, e.g., claim size distribution and ruin probability, 
age distribution and expected number of particles alive, service time 
distribution and stationary waiting time distribution, or innovation 
distribution and tail area for weighted averages. Mainly, these anal¬ 
yses are first order asymptotic results. Some second order results 
have been obtained. For compound or subordinated distributions, 
that is distributions of the form with F subexponen¬ 

tial, second order results have been obtained in Omey and Willekens 
(1987). For finite order convolutions, Baltrunas and Omey (1998) 
prove a second order result. 

Throughout this paper, we consider an infinite weighted sum 
where the A^’s are independent and identically dis¬ 
tributed random variables. We write F the distribution function 
of the XiS and G that of the weighted sum. 

Davis and Resnick (1988) consider the tail of G when F is in the 
domain of max-attraction of the Gumbel distribution, and belongs 
to a class slightly larger than the subexponential one. To state 
in a compact form the part of their result which we will use, and 
because we will be mostly interested in distributions and not in 
random variables, we introduce the following notation. If X has 
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distribution F, we write M^F the distribution of cX. Thus, whenever 
c is positive, McF = F{-/c), while MqF is the distribution function 
of the point mass at 0. Since it acts on functions, we call Me a 
multiplication operator. The distribution of Xliez is the infinite 
weighted convolution -ki^zMc^F. 

In this introduction, we also write C(i) the largest weight, and 
its multiplicity in the sequence of weights (ci). 

Theorem. (Davis and Resnick, 1988, Proposition 1.3). Let F he 
a subexponential distribution supported on the nonnegative half-line 
and in the domain of max-attraction of the Gumbel distribution. 
Assume furthermore that the Ci ’s are nonnegative and that 

cf is finite for some positive 6 less than 1. (1-1) 


Then 

( 1 . 2 ) 

We remark that Geluk and De Vries (2004) show a more general 
version of the above result; see also Chen, Ng, Tang (2005). 

As the above result shows, in this case the largest weight con¬ 
trols the first order asymptotic equivalence. This stands in contrast 
to what happens in the regular variation heavy tail case, where all 
weights are on an equal footing and appear in the first order asymp¬ 
totic equivalence; see Resnick (1987, p.227). One might propose that 
the next term in an expansion would be, say, k 2 Mc^.^^F where C( 2 ) is 
the second largest coefficient and ^2 equals the multiplicity of C( 2 ) 
in the sequence (c^). This intuition is correct for a narrow range of 
distributions. In section 2, we consider three regimes within the light 
subexponential tail one. In the heaviest of these regimes, the above 
intuition provides the right answer for the expansion; see Theorem 
2.3.2. In the other two cases, the missing element in our discussion 
makes its appearance, namely, derivatives of the underlying distri¬ 
bution. In fact, our formalism will not only provide the second order 
term, but, more generally, arbitrary many terms provided the distri¬ 
bution F is smooth enough. 

We conclude this introduction by noting sufficient conditions 
for a distribution to be subexponential have been obtained by 
several authors. Simple sufficient conditions in terms of the hazard 
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function — logF are given in Pitman (1980), Teugels (1975) and 
Kliippelberg (1988). Another sufficient condition on F is given by 
Goldie (1978). Finally, Goldie and Kliippelberg (1998) and Goldie 
and Resnick (1988) give sufficient conditions for a distribution to 
be both subexponential and in the domain of max-attraction of the 
Gumbel distribution. 

The paper is organized as follows. The next section contains our 
main results. A brief illustration of these results is given in section 
3, mostly to show that the main discussion in Barbe and McCormick 
(2005) carries over to the setting of the cnrrent paper. The proofs 
of the main resnlts are in the last section. 

2. Main results. In order to state our asymptotic expansions, we 
need first to discnss more precisely the class of distribution fnnctions 
which we will nse, and some algebraic tools which will allow us to 
both state and compnte these expansions. This will be done in the 
next two sub-sections. Onr main resnlts are presented in the third 
subsection. 

Results for some specific random weights require extra develop¬ 
ments and will be reported elsewhere. 

2.1. Hazard rate and smoothly varying functions. Consider 
for the time being a distribntion function F supported on the 
nonnegative half-line. Assuming it exists nltimately, the hazard rate 
h = F'/F will be of primary importance. Our analysis will show that 
how the hazard rate fnnction compares with logt asymptotically 
determines the form of the asymptotic expansion for tail areas. 

Writing the tail fnnction as 

F(t) = F(to) exp^—y h{u) du^ , 

we see that if h{t) ^ a/t then F is regnlarly varying with index a. 
On the other hand, if limt^oo h(t) = a, then F{t) = 
has a behavior close to that of an exponential distribution. Thus, 
in order to stndy tail area behavior for weighted convolutions in the 
lighter than regularly varying case bnt heavier than exponential case, 
we are led to consider hazard rates which satisfy 

h is regnlarly varying, 

(. 1 . 1 ) 

lim^^oo th{t) = -|-oo and limt^oo h{t) = 0 . 


4 



All the distribution functions satisfying assumption (2.1.1) are 
rapidly varying; indeed, writing Id for the identity function on the 
real line, for any M, the hazard rate is ultimately more than M/Id. 
Therefore, for t large enough and positive a less than 1, we have 

= expT— [ h{u) du) ^ expf— [ M/udu) = , 

F{at) ^ Jat ' ^ Jat ' 

showing that lim^^oo F(f)/F(at) = 0. It follows that if 0 < a < 6 
then MaF = o{MbF). 

The representation of distribution functions in terms of their haz¬ 
ard rate is closely related to the representation for the class T of 
De Haan (1970); see Bingham, Goldie and Teugels (1989, §3.10). In 
particular, if (2.1.1) holds, then F{t + x/h{t)^/F{t) = 

and F belongs to the domain of max-attraction of the Gumbel dis¬ 
tribution (see Bingham, Goldie, Teugels, 1989, Theorem 8.13.4). 
Moreover, if the index of regular variation of the hazard rate is 
negative, then Pitman’s (1980) criterion implies that F is subex¬ 
ponential. Consequently, if the hazard rate satisfies (2.1.1) and has 
a negative index of regular variation, then the corresponding distri¬ 
bution function satisfies the assumptions of Proposition 1.3 in Davis 
and Resnick (1988), and (1.2) holds provided (1.1) holds as well. 

For most linear time series models considered in practice, the 
sequence has an infinite order moving average representation (see 
Brockwell and Davis, 1987, Theorem 3.1.3). In that setting, F is the 
innovation distribution, and it is often unrealistic to assume that it is 
supported on the nonnegative half-line. A more common assumption 
then, is that of tail balance, asserting that F/M^iF has a positive 
finite limit at infinity. This assumption is not necessary for our 
results to hold, and some possible variations should be clear; but 
it makes results easier to state, and so we will use it. We then 
write h- the hazard rate of /W_iF when it exists. We then have the 
representation for the lower tail of F given by 

F{—t) = F{—to)exp(^—J h-{u)duj, 

for t at least to- Note that if the distribution is tail balanced, then 
J^^(h — h-){u) du has a limit at infinity. This implies that if both h 
and /i_ are regularly varying with index greater than —1, then they 
are asymptotically equivalent. 
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Throughout this paper we adopt, without further mention, the 
following convention. 

Convention. If a distribution function F is tail balanced and 
its hazard rate h satistifies a property, then h- satisfies the same 
property. To emphazise this convention, we say that F is strongly 
tail balanced. 

For instance, if the property is ’/i has a continuous m-th deriva¬ 
tive’, we mean that h- has also a continuous m-th derivative. 

Assuming that the hazard rate is regularly varying is not always 
sufficient to obtain higher order results, in particular because higher 
order derivatives of the hazard rate will be involved. This moti¬ 
vates the following definition, introduced in Barbe and McCormick 
(2005). It somewhat extends that of a smoothly varying function as 
presented in Bingham, Goldie and Teugels (1989). 

Definition. A function h is smoothly varying of index a and order 
m if it is ultimately m-times continuously differentiable and the m-th 
derivative is regularly varying of index a — m. We write SRa^m 
the set of all such functions. 

Recall that Id is the identity function on the real line. For any 
real number a and any nonnegative integer m, we write {a)m for 
a{a — 1)... (a — m -|- 1). Note that if h is smoothly varying of index 
a and order m, then for any nonnegative integer k at most m, the 
asymptotic equivalence h^^^ ~ (a)fc/i/Id^ holds. 

2.2. Laplace characters. To effectively calculate asymptotic 
expansions, we need to recall the algebraic method introduced in 
Barbe and McCormick (2005) in their analysis of tail areas for linear 
and related processes based on innovations with regularly varying 
tails. We recall the terminology and concepts from that paper for 
which we shall have use. We write D the derivation operator, that 
is, if is a differentiable function. Tig is its derivative. We write D® 
the identity operator, and D* = for any positive integer i. 

For a distribution function F, we write fiF,i its i-th moment. 

Because they will appear in our asymptotic expansions, the basic 
algebraic objects that we need are the Laplace characters. 
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Definition. (Barbe and McCormick, 2005). Let F be a distribu¬ 
tion function having at least m moments. Its Laplace character of 
order m is the differential operator 


-F,m = 

O^i^m 


iF.m. — 


{-ly 


i\ 


-hF, 




The algebraic role of these operators derives from their property as 
morphisms from the convolution semi-group of distribution functions 
with at least m moments to the space of differential operators modulo 
the ideal generated by More specifically define a bilinear 

product on (equivalence classes of) differential operators modulo the 
ideal generated by denoted o, by 

D* o D-^' = I ifi + j^m, 

\ 0 otherwise. 

This is just an explicit expression for composition of differential 
operators modulo the ideal generated by Then, it is a simple 

matter, involving only the definition of convolution and the binomial 
formula, to show the useful relation 

We refer to Barbe and McCormick (2005) for a more extensive 
discussion of Laplace characters and their practical use. 

2.3. Asymptotic expansions. To write concise formulas, we 
need to introduce a bit more specialized notation, tj. This symbol is 
used to delete a factor from a convolution product. More precisely, 
if LI = Ki -k K 2 is a convolution product of distributions Ki and 
K 2 , then H\\K 2 = Ki. In particular G'^Mc^F is the distribution of 

J 2 jez\{i} ^3^3■ 

To allow weights and random variables to be negative and still 
write concise statements, we introduce a partial ordering on real 
numbers, reflecting some form of tail behavior. If a and b are two 
real numbers, we write a ^ b if MaF = o{Mi,F). If there exists a 
positive e such that cMaF if MbF ^ e ^ MaF ultimately, then we say 
that a is equivalent to b and write a = b; provided the denominator 
does not vanish ultimately, this means that the ratio MaF/ MbF is 
ultimately bounded away from 0 and infinity. We also write a Fb if 
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and only if a ^ 6 or a = 6 . Naturally, we define a y b as b ~< a and 
a y b as b y a. 

Since F is rapidly varying, if both a and b are positive, then a y b 
if a < 6 , while a = 6 if and only if a = 6 . When a and b have 
different signs, the ordering introduced depends on how both tails of 
F compare. In particular, if F vanishes in a neighborhood of —oo, 
then a = 0 whenever a is negative, and ^ defines a total ordering 
on the real numbers. Note that the ordering is also total when F is 
tail-balanced or vanishes in a neighborhood of —oo. 

We use the notation to represent the sequence 

in nonincreasing ^ order and without multiplicity; by the latter, 
we mean that C(i) is strictly decreasing for the order, except in 
the case that it eventually reaches the minimal element 0 when it 
stays identically constant 0 from that point onward. For instance, if 
Cl = C 2 = 1 while cs = 1/2 and 04 = 05 = 1/3 and all the other Ci’s 
are nonnegative less than 1/3, then C(i) = 1, while C( 2 ) = 1/2 and 
0 ( 3 ) = 1/3 and so on. Note that any sequence with limit 0 

at infinity can be put in nonincreasing order indexed by the positive 
integers. We say that an element Ci of the sequence is maximal if it 
is equivalent to C(i). 

Since we will not restrict the sequence {ci)i^z to be nonnegative, 
we need to replace condition ( 1 . 1 ) by 

Icil"^ is finite for some positive S less than 1 . (2.3.1) 


The following convention will be convenient. 

Convention. We say that the ’standard conditions’ hold if (2.3.1) 
holds and either F vanishes on a neighborhood of —00 or is strongly 
tail balanced. 

Minor changes in our proof show that the results of this paper 
remain true if instead of requiring F to vanish in a neighborhood 
of —00 we impose that dF(x} is finite whenever we use the 

Laplace character of order m of T in these results. 

Our first result covers most light subexponential distributions 
used in applications, that is those for which — logF is regularly 
varying of index positive and less than 1. In particular, it includes 
Weibull distributions; however, the log-normal distribution is ex¬ 
cluded. 



Theorem 2.3.1. Assume that the standard conditions hold and 
(i) h{t) > t~^ logt, 

(a) h is smoothly varying of negative index and order m, 

Then 

G= + oih^M:^). 

To convey the meaning of this formula, let us consider the case 
when the sequence Ci is nonnegative, indexed by the positive integers 
— that is Ci vanishes if i is negative — decreasing, with ci = 1. 
Writing H = G'^F, that is the distribution of X]i >2 ^i^ii formula 
in Theorem 2.3.1 asserts that for any k at most m, 

G = F- iiH,iF' + ^F" + • • • + . 

2 k\ 

More generally, for a sequence Ci of arbitrary sign with maximal 
element C(i) = 1, and assuming that F is tail-balanced, Theorem 
2.3.1 yields the following. Let (respectively v-) be the number 
of maximal positive (respectively negative) cfs which are maximal 
for the ^ ordering. Since F and /W_iF are rapidly varying and 
asymptotically of the same order, the maximal positive c^’s, if they 
exist, are all equal, and equal to minus the maximal negative cfs if 
they exist. Consequently, Theorem 2.3.1 asserts that for k at most 
m, 


Qualitatively, this result is very much in agreement with Davis 
and Resnick’s result, showing that the largest weight drives the 
asymptotic behavior of G. 

For practical purposes. Theorem 2.3.1 is applicable over a broad 
range of Weibull-like subexponential distributions. However, its 
scope does not include the useful subexponential distribution given 
by the log-normal distribution. Moreover, Theorem 2.3.1 stands in 
contrast with the situation where F is regularly varying, where each 
weight contributes to a term of the asymptotic expansion (Barbe and 
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McCormick, 2005). The purpose of the remainder of this section is 
to show the delicate transition between regularly varying tails and 
those covered by Theorem 2.3.1. 

Theorem 2.3.1 contains the assumption that the hazard rate is 
asymptotically much larger than Id~^ log, excluding distributions 
with tail when a is greater than 1 and at most 2. The next 

result covers this class of distributions with a less than 2 (note 2 is 
not allowed yet!). It shows that this narrow class of distributions, as 
far as the asymptotic behavior of weighted convolution goes, retains 
some features of regularly varying tails, yet has rapidly varying tails. 

Theorem 2.3.2. Assume that the standard conditions and (2.1.1) 
hold, and 

(i) h{t) = o{t~^logt) as t tends to infinity 

(a) th{t)^/h[l/h{t)) < oo, (2.3.2) 

Then, for any positive integer m, 

G= 

Z.Cj 


Setting r]{t) = th{t), condition (2.3.2) asserts that the supremum 
limit of rj{t)/ri{l/h{t)^ is finite. This is the case if rj is of order log“. 

Note that under the assumptions of Theorem 2.3.2, the index of 
regular variation of the hazard rate is —1. In the conclusion of this 
theorem, the largest term of the expansion is driven by the largest 
coefficient c^, the second order term is driven by the second largest 
coefficient, and so on. Although formally the expansion has a similar 
form as a first order equivalence in the case of distributions with 
regularly varying tails, namely Yhie-L^ciF, here, because the tails 
are rapidly varying, nonequivalent weights produce asymptotically 
distinct contributions to the expansion. When F is rapidly varying, 
the MciF are ordered by tail dominance, and this yields formally the 
formula presented in Theorem 2.3.2. 

Our last result yields some very unexpected behavior and fills the 
gap between Theorems 2.3.1 and 2.3.2. It covers the intermediate 
range where the hazard rate is exactly of order t~^ log t and therefore 
logT is exactly of order — log^. In particular, it covers a distribution 
often used in practice: the log-normal. A close look at the conclusion 
of that theorem (particularly examining limiting cases as A tends to 0 
or infinity) reveals a natural progression in the terms involved in the 
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expansions in Theorems 2.3.1 and 2.3.2 with only terms equivalent 
to C(i) in the former result and as many as desired in the latter result 
(viz. A = oo and A = 0). 

In this result [-J denotes the integer part, that is the largest 
integer at most equal to the argument. 

Theorem 2.3.3. Assume that the standard conditions hold and 
(i) h{t) Xt ^ log t for some positive X, 

(a) h is smoothly varying (necessarily of index and order m. 
Then, for any positive integer k at most m, 

G = ^ C(i)e ^ }^Gtl/Wc^F,fc+[A log(|ci/c(i) |)J 

iez 

+ oih’^Mf^) . 

With regard to this theorem, we note that the presence of the integer 
part serves to indicate that for small Ci, the order of the Laplace 
character may be reduced. The indicator function ensures that no 
negative order appears as an order to a Laplace character. Since 
differentiation and multiplication operators change the asymptotic 
order of their argument in different ways, it may happen, for 
example, that is of higher order than yielding a 

case where the second largest coefficient has dominance over the 
largest coefficient in determining the second order term. This will 
be illustrated in the next section. 

3. Examples. The purpose of this section is to illustrate briefly our 
three results, and show that, as was done more lengthily in Barbe 
and McCormick (2005) for regularly varying tails, the asymptotic 
scale in which an expansion is written influences which terms of our 
results should be kept. 

The following four examples illustrate successively Theorems 
2.3.2, 2.3.3 and 2.3.1. Throughout these examples, we consider a 
sequence of weights such that Ci is 0 if i is nonpositive and the se¬ 
quence Ci is nonnegative and strictly decreasing for i positive. We 
also take ci to be 1. Finally, the summability condition (2.3.1) is 
assumed to hold. In the hfth example, we derive a three terms ex¬ 
pansion for the marginal distribution of some linear processes. 

Example 1. Dehne the functions 
ei(t) = exp(—log^^^ f) and e 2 (t) = ei(t) exp(—log^^^ t), 
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so that 62 = o(ei) at infinity. Consider a distribution function F 
such that F = 61 + 62 ultimately. Since + ~ 61 , Davis and Resnick 
(1988) shows that G ^ F ~ ci, because our sequence Ci is decreasing 
with Cl = 1. What is the second order? 

Note that for any positive number a less than 1, 

logei(t/a) - log 62 (t) ~ -(logt)^/^ logo. 

Consequently, MaCi = 0 ( 62 ). Hence, for any i at least 2 , we have 
Mc-F = 0 ( 62 ). It follows from Theorem 2.3.2 that the relation 
F = 61 + 62 , holding ultimately, shows that F provides not just 
a first order equivalence to G but actually a second order result in 
that G = Cl + 62 + 0 ( 62 ). 

Example 2. Let 0 be a positive real number. Take hit') = 20t“^ logt, 
so that F = exp(—01og^) on [l,oo). Define 7 = X]i 7 i+- Note 
7 is nonnegative and finite by our assumptions. We see that 
= 7 +f,i- Since ci = 1, Theorem 2.3.3 shows that a two 


terms expansion is either 



case 1 : 

F - 7 /Uf,iF' 

if 

II 

case 2 : 

F + M,,F 

if 

/We,F»F, 

case 3 : 

F + Mc^F - 7 +f,iF 

if 

X F' 


We are in the first case if 62 ^ and in case 2 otherwise. Case 

3 cannot occur for this specific distribution. But case 3 may occur 
if, for instance, F{t) = exp(—log^ t + loglogt) ultimately! Note 
also that for the log-normal, h{t) = logt +corresponding to 
9 = 1 / 2 . 

Example 3. Assume again that ci = 1 and moreover that ki = 1. We 
write p,i for We suppose that the assumptions of Theorem 

2.3.1 are satisfied for m = 3. The conclusion of Theorem 2.3.1 yields 
a seemingly four terms expansion, 

G = F- /riF + ^F" - ^F'" + o(FF) . (3.0.1) 

2 6 

However, this formula hides the true number of significant terms in 
the expansion. To clarify this remark, we further investigate the 
meaning of (3.0.1) and calculate 

F' = -hF, 

F" = i-h' + h^)F, 

F'" = i-h" + 2>h'h-h^)F. 
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It is then natural to seek an expansion in the scale containing all the 
functions involved in these derivatives, that is 

F, hF, h^F, h^F, h'F, hh'F, h"F 

which are not o{h^F). Then, for example, a three terms expansion 
in that scale can be obtained from (3.0.1)and is given by 

G = F + nihF + + o{h^F). 

In similar fashion, from (3.0.1), using the expression for the deriva¬ 
tives of F, we calculate that a four terms expansion, where now all 
terms in the formulas below are significant, is either 

G = FF nihF + ^h^F - ^h'F + o{h'F) 

or 

G = F + fiihF + ^h^F - ^h^F + oih^F ), 

2 o 

or 


G = F + nihF + ^h^F - ^h'F - ^h^F + o{h^F ), 

2 2 6 

the first one occurring if = o{h'), the second one occurring if 
h' = o{h^) and the third one occurring if both and h' are of the 
same asymptotic order. Note that if the index of regular variation 
of — logF is less than 1/2 then we are in the first case, while if it is 
more than 1/2 we are in the second case. When this index is equal 
to 1 /2 we can be in any of the three cases a priori. 

Example /. The purpose of this example is to show that some 
cancellation may occur, analogous to that observed for regularly 
varying tails in Barbe and McCormick (2005). 

Consider a distribution for which ultimately 

F{t) = exp(—log^^^ t) — exp(—log^^^(2t)) . 

Since 

log^/^ ( 2 t) = log^/^ t + ^ log^/^ t log 2{l + o(l)) , 
we see that 

F{t) ~ exp(— log^/^ t). 
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Assume that Ci = 1 and C 2 = 1/2. Theorem 2.3.2 with m = 2 yields 
G{t) = F{t) + F{2t) + o(F{2t)) , 

apparently giving a two terms expansion. However, note the cancel¬ 
lation between the second order term of F{t) and the leading term of 
F{2t), so that the above expansion gives only the one term expansion 

G{t) = exp(- log^/^ t) + o[F{2t)) . 

It is clear from this example that a construction similar to that 
done at the end of section 3.2 in Barbe and McCormick (2005) can 
be made: regardless whether we apply any of the three theorems of 
the previous section, for any fixed m, we can find a distribution and 
a sequence of weights such that the asymptotic expansion of G in the 
theorem gives in fact a one term expansion, because of cancellations. 
Hence, we see once more that obtaining a second order term for an 
infinite convolution is not that of obtaining a second order expansion, 
but that of obtaining arbitrarily accurate expansions. 

Example 5. Consider a stationary linear process with innovations 
having a symmetric distribution F. The stationary distribution, 
G, is that of some proper constants c^. Since F is 

symmetric, the odd moments of G vanish. 

We assume that the sequence (c^) has a unique maximal element 
in the ^ ordering, equal to 1 say. 

To calculate the first two even moments, for any positive integer 
k, write Cl = cf • Then 

= MGtlF,2 = CIpf,2 

l‘Gtl/Wc(jjF,4 = MGtlF,4 = 3(C2^ — CI)Pf2 + CIpf,A ■ 

If F obeys the assumptions of Theorem 2.3.1 with m = 6 say, we 
have 

G = F+ jPP' + 3(C’2-^ - QK.a + C4>f.« p(4) + o,y(6)) 

This formula applies for instance to symmetric Weibull type densities 
Q,g-kr/2r(l/a;). It applies also to densities which are ultimately 
equal to x^e~^ ; compare with Rootzen (1986) who obtained first 
order results. 
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4. Proofs. The proofs are based on the writing of the convolutions 
in terms of operators modelled after those introduced in Barbe and 
McCormick (2005). For every distribution function K, define the 
operator 

ptl2 

TKf{t)= fit-x)dKix). 

J —OO 

Then, 

Hi.K = TkH + ThK + M2(HK) . (4.0.1) 

Throughout the proofs, we write F) for Mc^F. Without loss of 
generality, we assume that the index set of the sequence (c^) is N 
and not Z. 

Define Gn = Fi -k ■ ■ ■ -k F^. By induction, formula (4.0.1) yields 



* —fc+2 —fc + l 

2^k^n 


+ M2{Gn-l Fn) + T . 

•• Tp'^_j^_^^M2{Gn — kFn — k-\-l 

2^k^n 



(4.0.2) 


Note. Throughout this section, we assume that (2.1.1) holds, 
possibly without mentioning it. 

4.1. Some lemmas. Preliminary lemmas are presented under not 
necessarily sharpest conditions. Use of these lemmas appears mainly 
in the proof of Theorem 2.3.1. Our first lemma gives the exact order 
of derivatives of F. 


Lemma 4.1.1. Assume that h is in SRa-i,m for some a nonneg¬ 
ative and less than 1. If a vanishes, assume further that th{t) tends 
to infinity at infinity. For any nonnegative integer k at most m, 


Proof. Using the Faa di Bruno formula (see Roman, 1980, p.809 for 
the formulation that we are using) equals 



l^i^k 


E 

niH- ^ni = k 


k\ 

nfi... .nil 
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Since h is smoothly varying of order m and k is at most m, for any 
index rij involved in the sum, ~ {a - l)nj-ih/l(i^^ \ As 

usual in the theory of regular variation, when (a — vanishes, 

such asymptotic equivalence must be read as = o(/i/Id”^“^). 

Thus, 

hin,-l) ^K-1) ^ _ l)^^_^/,Yld"-* , 

where we use ni + • • • + = k. The function /i*/Id^~*, that is 

(Id/i)*/Id^, is regularly varying of index ia — k. Among such terms, 
the one with the largest order is obtained when i = k, even if a 
vanishes, for in this case Idh tends to infinity at infinity. This forces 
all the Ui’s to be 1 and yields the result. ■ 

The next lemma is a Potter type bound which we have already 
stated and proved in section 2. 

Lemma 4.1.2. Let M be an arbitrary positive number. There 
exists ti such that for any t at least ti and any X at least 1, 

F{Xt)/F{t) ^ . 


Our next lemma shows that for t large enough, the function 
F{t — ■)F is nonincreasing on some interval [M, t/2], where M does 
not depend on t. 


Lemma 4.1.3. Assume that h is in SRa-i,m for some m at least 
1. There exists some positive ti and M such that for any t at least 
ti, the function x ^ F{t — x)F(x) is nonincreasing in [M,t/2]. 


Proof. Let (5 be a positive real number such that a —1+(5 is negative. 
Lemma 2.2.4 in Barbe and McCormick (2005) shows that there exist 
ti and M such that for any t at least ti and x in [M,t/2], 


/ /t \ 


This upper bound is nonpositive in the given range of x, implying 
the result. ■ 
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The following result will be instrumental. 


Lemma 4.1.4. Let K he a eumulative distribution function on the 
nonnegative half-line, such that for some M positive, K ^ MF. For 
any interval [a, 6] and any nonnegative function f continuous and 
nondecreasing on [a,b]. 



fdK^MU fdF + fFib)] 


[a, 6] 


We will sometimes use this lemma in its limiting form, as b tends 
to infinity. In particular, if both F and K are supported by the 
nonnegative half-line, it implies that HK,i ^ for any positive 

integer i. 


Proof. Assume that / is differentiable. An integration by parts 
yields 



fdK = fK{a)-fK{b) + 



f K{x) dx. 


Hence, the left hand side of (4.1.1) is at most M times 


(4.1.1) 


fF{a)+ [ f'F{x)dx. 

J [a,b] 

But taking AT to be F in (4.1.1) yields that this sum is at most 



fdF + fF{b). 


This prove the lemma when / is differentiable. When it is only 
continuous, use that the differentiable functions are dense in the 
continuous ones in (L^[a, 6 ], dK). ■ 


Combining the two previous lemmas yields the following asymptotic 
rate of decay on the tail of moment integrals for a class of distribu¬ 
tions. 


Lemma 4.1.5. Let K be a distribution function such that K = 
0{F). Assume that the hazard rate function h of F satisfies (2.1.1). 
For any nonnegative integers j and n. 



dK{x) 


o{t-n • 


17 



Proof. There exists some positive e such that eK ^ F ultimately. 
Using Lemma 4.1.4, we have for t large enough 


x^dK{x)^ / x^ dF{x) 


= FF{t)+jF I ^F{ts)ds 


^FF{t)+jFF{t) I s-^ds. 


This bound is rapidly varying, hence it is o{t "■) for any n. 


4.2. Proof of Theorems 2.3.1 and 2.3.3. In this subsection, we 
first prove Theorems 2.3.1 and 2.3.3 with the extra assumption that 
the weights c^’s are nonnegative. We will remove this assumption 
afterwards. So, from now on, until stated otherwise, we assume that 
the weights are nonnegative. 

We start by proving a lemma which allows us to neglect terms 
arising from the multiplication operator in formula (4.0.2). 


Lemma 4.2.1. Assume that the hazard rate is regularly varying 
of negative index and that liminf^^oo th{t)/\ogt is positive (possibly 
infinite). Then, for any k nonnegative, M 2 F =o{h^F). 


Proof. Using the representation for F, for any t at least to, 


/i^F 


(t) = F(to) exp 


nt/2 

/ h{u) du+ h{u) du — k log 
to t j‘2 


Let e be a positive real number. Using that h is regularly varying of 
index a — 1 for some a less than 1, we have ultimately 




du^ (1 + 0 ( 1 )) 


The right hand side above is less than some —ath{t) for some positive 
a by taking e small enough. 

If a is positive, and less than 1, then the result is clear since 
logh{t) = o{th{t)). 
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If a vanishes, a can be taken as large as desired by taking e small 
enough. By assumption th{t) is ultimately at least like some positive 
multiple of log t, implying that log h{t) grows at most as a multiple of 
logt; moreover, th{t) is is ultimately at least some positive multiple 
of logt. This proves the lemma. ■ 

For any positive real number M, we write B{F, M) the set of all 
distribution functions K such that both K and M^iK are less than 
MF on the nonnegative half-line. 

Our next lemma shows that remainder terms stay negligible under 
the action of some operators. 


Lemma 4.2.2. Assume that h is regularly varying with negative 
index of regular variation and F is subexponential. Let K be a 
distribution function in some B{F,M). If g is any function which 
is o{h"^F) then, 

TKg = o{h^F ). 

Proof. Let e be a positive real number. Let t be large enough so 
that ^ eh^F on [t/2,oo). Theorem 1.5.3 in Bingham, Goldie 
and Teugels (1987) implies that h{t) is asymptotically equivalent to 
suPs^t h{s). For t large enough, we obtain for some positive p 


.t/2 _ 

\TKg{t)\ife j {h'^F){t - x)dK{x) 

J — oo 

1-112 _ 

ffe2Ph^{t) / F{t-x)dK{x). 


f /2 

Lemma 4.1.4 implies that Jq' F(t —x) diL(x) is less than a constant 
times TpF{t) -|- M 2 F^{t), which is at most a constant times F*^{t), 
while J^^F{t — x)dK{x) is at most F{t). The conclusion follows 
using subexponentiality of F. ■ 

The next lemma is at the heart of both Theorems 2.3.1 and 
2.3.3. We will use it to approximate of some T operators by Laplace 
characters, when acting on derivatives of F. In the proof, we will 
use the absolute moments 

\h\K,i = J kr diL(x). 
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Lemma 4.2.3. For any fixed integer p at most m and any positive 
M, 


lim lim sup sup 

^ ^0 t —^oo 


h^F{t) 


LK,m-pF(-P\t) 


= 0 , 


where the supremum is taken over all distribution funetions K 
belonging to B{F,M). 


Proof. Taylor’s formula with remainder term asserts that 


p —1 *^’ 

px rn-p-l _ 

+ (-ir“M --F(-)(t-x + y)dy. 

Jo [m-p-iy. 

Integrating this equality between —S/h{t) and 6/h{t) with respect 
to dK{x) shows that 


fS/h(t) 


'-5/h(t) 




is at most 


E 


\F(P+F(t)\ 


\xj dK{x) 


O^j^m-p 

r-5/h(t) 

+ / 

J —S/h(t) 


—5/h{t),S/h{t)]^ 

(F(m) ^ + y) _ pim) ^ 

[m — p — 1)\ 


ym-p-l 


Combining Lemma 4.1.1 and the proof of 4.1.5, we see that for any 
positive (5, 


lim sup 

t^oo xeB{F,M) 


i[_s/h{t),5/h{t)Y 1 ^ 1 '^ AK{x) 
h"^F{t) 


= 0 . 


Let e be a positive number. For any fixed 6, Lemma 4.1.1 and 
regular variation of h show that uniformly in z in [—1,1] say, 
(t — 5z/h{t)) /F^'^\t) tends to e^^. Therefore, we can hnd <5 
such that, ultimately, for any z in [—1,1 ], 


1 -e^ 


F^^^ {t-Sz/h{t)) 


^ 1 + e . 
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Consequently, ultimately and uniformly over B(F,M) 


pS/hit) 

l-5/h{t) 


y 


m p 1 / ^(m) 


Jo {m-p-1)1 

pS/h{t) 

J-S/h{t)'J0 

^ (-lylK,m—p 


y 


m—p — 1 


— 1 ^ dy dK{x) 


{m — p — 1)! 


dy 


dK{x) 


Since e is arbitrary, Lemma 4.2.3 follows from Lemmas 4.1.1 and 
4.1.4. ■ 


As announced, the preceding lemma yields an approximation of 
some T operator by Laplace characters, which we now state and 
prove. 


Lemma 4.2.4. For any p at most m, 


lim sup 
KeB{F,M) 


{Tk - LK,m)F(-P^ 


{t) = 0. 


Proof. Given Lemma 4.2.3, it suffices to show that for any positive 
e and any S small enough. 



-x)dA(x) ^eh^F{t) 


(4.2.1) 


and 


r--S/h(t) 


F^P\t-x)dK{x) 


^ eh^^Fit) 


(4.2.2) 


ultimately and uniformly over B{F,M). Using Lemma 4.1.1 and 
convergence of the hazard rate to 0 at infinity, the left hand side of 
(4.2.1) is ultimately bounded by F{t — x) dK{x). By Lemma 

4.1.4, this is less than 


/■*/2 _ _ 

M / F{t - x) dF{x) + MF^{t/2). 

For t large enough, we substitute dF(x) by F'{x) dx in this bound. 
Thus, Applying Lemmas 4.1.1, 4.1.3 and 4.2.1, this upper bound is 
ultimately at most 


MtF{t - 5/h{t))F{5/h{t))+o{F^F{t)) . 
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This implies (4.2.1) because on one hand, F{t — 5/h{t)) F{t)e-\ 

and, on the other hand, F being rapidly varying, F(5/K) = o(/i^) 
for any positive p, and, in particular, IdF((5//i) = o{h'^). 

To prove (4.2.2), we replace F^p"> by h^F using Lemma 4.1.1. 
Since x is negative in the range of integration, we bound h'^{t — x) 
by supg>j/i™'(s), which in turn is equivalent to and bound 

F{t — x) by F{t). Then (4.2.2) follows, and this completes the proof 
of Lemma 4.2.4. ■ 


Note that in Lemma 4.2.4, we could replace F by MaF for any 
positive a, since the hazard rate of MaF is a~^h{a-) and has the 
same properties as h as far as the proof of the lemma is concerned. 

We now complete the proofs of both Theorems 2.3.3 and 2.3.1, 
assuming positivity of the weights. Without any loss of generality, 
we assume that C(i) is 1. Note that in Theorem 2.3.3 the hazard 
rate is regularly varying of index —1. While, in Theorem 2.3.1, the 
hazard rate is regularly varying of index a — 1 with 0 ^ a < 1. 

Recall formula (4.0.2). Again, using Davis and Resnick’s (1988) 
Proposition 1.3, Gn-kFn-k+i = 0{F'^). Thus, Lemma 4.2.1 implies 

M2(Ga-kFa-k + l) = 0{h^F) , 

and by induction, using Lemma 4.2.2, 

M 2 (Gn-iF^) + TF„^,^,M 2 (Ga-kFn-k+i) = o{h^F ). 

2^k^n 

It follows that under either the assumptions of Theorem 2.3.3 or 
Theorem 2.3.1, 

Gr.= TG„_,Fa+ Y TF„...TF„_,^,TG^^,F^-k + l+o{h^F). 
2^k^n 

Again, combining Davis and Resnick’s (1988) Proposition 1.3 with 
Lemma 4.2.4, we see that 

^Gn-k^n-k + l = l-G^_f.,mFn-k + l + o{h F) . 

Any Laplace character applied to some Fi yields a linear combination 
of derivatives of Fi. Therefore, using Lemmas 4.2.2 and 4.2.4 
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inductively, we see that 

Gn = y ^ ^Fn,m o ■ ■ ■ O ,m-^n-fc + l + o{h F) 

l^k^n 

= E l-G„^F^.,^,,mFr.-k+i + o{h-^F). (4.2.3) 

l<k<n 


Note in the first sum above the term indexed by A: = 1 is to be read 
as LG„_^,mFn- Write Fin for the distribution function of -ki^n^aF. 
We claim that there exists n large enough so that Hn = o{hF'F). 
Since our sequence (c^) is nonincreasing, Davis and Resnick’s (1988) 
Proposition 1.3 implies that Hn ^ Fn+i and 

(^) = ^^p(-/ /i(u)du-mlog/r(t)j 

/ c~“ — 1 \ 

= eyip\^—th{t) -(l + o(l)) — mlogh{t)j .(4.2.4) 

Thus, our claim is now evident. Further, (4.2.4) shows that under 
the assumptions of Theorem 2.3.1, in order that Hn = o{h'^F), it 
suffices that n be chosen large enough so that c„ +1 < 1 . 

Next, note that G = Gn * Hn, so that 

G = Tfj^Gn + Tc^Hn + M2{Hn Gn) • 

Since Hn = o{h^F) and Davis and Resnick’s (1988) Proposition 1.3 
asserts that Gn ^ F, Lemmas 4.2.2 and 4.2.1 show that 

G = Tn^Gn + o{h^F ). 

Thus, using Lemmas 4.2.4, 4.2.2 and formula (4.2.3), we obtain 
G= E LGm,mFk + o{h^F). 

l^k^n 

This proves Theorem 2.3.1 since we can take n to be the first 
integer for which c^+i is less than 1 , that is to be ki. 

This proves Theorem 2.3.3 as well for the following reason. First, 
recall that under the assumptions of that theorem, a vanishes and 
in (4.2.4) we need to read (c““i — l)/a as equal to — logc„+i. Then, 
since h{t) ~ Af^logt, equality (4.2.4) shows that in order to have 
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Fn+i = o{h"^F), we must have (—Alog l/c„_l_i + m) negative, that 
is Cn+i < Hence we have proved that 

G = I{ ^ + oih^F). 

Finally, since Fn'’ is of order FF^, it is of order smaller than 
K^F as soon as Cn < that is i > m + Alogc„. Thus, 

the differential operators in the Laplace characters only make a 
nonnegligible contribution for 0 ^ i ^ m + [A log c^J. Theorem 
2.3.3 follows. 

By an obvious symmetry, Theorems 2.3.1 and 2.3.3 follow if all 
the q’s are negative and F is strongly tail balanced. 

To remove the sign restriction on the constants, we write H 
(respectively K) for the distribution of 

CiXi (respectively CiXi ) . 

Then G is F[ -k K and 

G = ThK + TkH + M2(HK) . 

Now that we have obtained the asymptotic expansions for H and K, 
we conclude with nearly identical arguments. 

4.3. Proof of Theorem 2.3.2. We begin by proving a lemma 
which allows us to neglect terms arising from the multiplication 
operator in formula (4.0.1). Recall that under the assumptions of 
Theorem 2.3.2, the hazard rate is regularly varying of index —1. 

Lemma 4.3.1. Assume that h is regularly vaying of index —1 and 
(2.1.1) holds. For any positive numbers a, b, and c, 


MaF MhF = o{M,F). 


Proof. Since x i—> M^F is nondecreasing on (0, oo), it suffices to 
prove the lemma when 0 < c ^ a ^ 6. Using the representation of 
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F in terms of its hazard rate, we have for any positive 6 less than 1 
and for any t at least to, 


/M F MhF \ _ 

log( ^ / H'^)du- / h(u)du + logF(fo) 

^ Met / Jt/a J&t/b 

= th{t){\og{a/c) + log(5)(l + o(l)) . 

Taking 5 small enough ensures that \og{5a/c) is negative. The result 
follows since th{t) tends to infinity with t. ■ 

Our next lemma is the essential step in the proof. It shows that 
some T operators are close to the identity. 

Lemma 4.3.2. Let M be a positive number, and let K he a 
distribution function in B{F,M). Then, under the assumptions of 
Theorem 2.3.2, for any positive real numbers a and b, 

{Tk = oilVkF) 


at infinity. 

Proof. Let (5 be a positive number. Consider first the integral 
.5/h{t) -F{t/a) 


I —/ , N ■ dK{x) 

Ls/m F{t/b) 

_F{t/a) /F((t-x)/a) 


F{t/b) J-s/h{t)^ F{t/a) 

If the absolute value of x is at most S/h{t), then 


-l)dit:(x). (4.3.1) 


F{{t-x)/a) 

F{t/a) 


= exp 


— I a 


' 11 a 


h{u) du 


/ Kui-^f h{t{l + v)/a) 

= exp(— h{t/a) / -—— -du 

V a Jo h{t/a) 

is ultimately between exp(—2/i(t/a)|x|/a) and exp(2/i(t/a)|x|/a) 
Therefore, if 5 is small enough, 

,F{(t-x)/a) \ 4 

(— FlTr\ - l) dK{x) ^-h{t/a) \x\dK{x) 

-6/h(t)^ Tit/a) / a J_s/h(t) 


—5/h{t) 

4 

^ -h{t/a)\n\K,i ■ 
a 
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It follows that the absolute value of the right hand side in (4.3.1) is 
ultimately at most 


AF{t/a) 

aF{t/b) 


h{t/a)\iJ.\K,i • 


(4.3.2) 


Recall that h is regularly varying. So, up to the 4:\fi\K,i/a factor, 
this upper bound is 


exp 


(^~ /i(u) du + logh(t) + 0(1)^ 

= exp^-th(t) log(6/a)(l + o(l)) + logh(t) + 0(1)^ 


Since h{t) = o{t~^logt) and th{t) tends to infinity, th{t) = 
o[logh{t)). So (4.3.2) tends to 0 and so does the left hand side 
in (4.3.1). 

It remains to show that 



/ f{ )dK{x) = o{F{t/h)) 

Js/h(t) V a / 

/ f( - )dK{x) = o{F{t/b)) 

J — oo \ CL ' 

as well as 

F{t/a)K{5/h{t)) = o{F{t/b)) , 

and 

F{t/a)K{-S/h{t)) = o{F{t/b)) . 


Because of our assumption that both K and M^iK are 0(F), and 
the function x i—> F(^{t — x)/a) is nondecreasing , these assertions 
are implied by 

F{t/2^F{5/h{t)) ^ F{l/h{t)) ^ 

F{t/b) F{t/b) h{t) 

This last estimate holds because (4.3.2) tends to 0 at infinity, F is 
rapidly varying and h is regularly varying. ■ 


From Proposition 1.3 in Davis and Resnick (1988) and using 
Lemma 4.3.1, we see that for any positive S, 

M2(Gn-kFr.-k+i) X M2(MZF = o{M^) . 
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But then, using the definition of the T operators, 


TFr,-.k +2 ^2{Gn-kFn-k + l) ^ ^2 ^2 (G'n-fc-f"n-fc + l) — o{l\/l 2 sF) , 
and by induction, taking S small enough, we see that 
/W2 (G„-i FJ + = o(F„). 

2^k^n 

Therefore, we have 

Gn=Tc„^:^Fn+ ^ TG„_fcF„-fc+i + o(F„). 

2^k^n 

We apply Lemma 4.3.2 to obtain 

FCn — kFri — k + l — Fn — k+1 T o(F^) . 

Then, another use of Lemma 4.3.2 shows that 
FFn-k+2^iFn) = o{Fn) , 

and proceeding by induction, we obtain the expansion of Gn, 

Gn = Fn-k+l + o{Fn) ■ 

To finish the proof, let Fin denote the distribution function -ki^nFi, 
so that G = Gn Hn- Then, representation (4.0.1) gives 

G = Tjj^Gn + Tc^Hn + M2Gn M2F[n ■ 


Up to increasing n we can assume that c^+i < Cn- Note that by 
Proposition 1.3 in Davis and Resnick (1988), Hn is of order F„+i, 
which is o{Fn)- Then, by Lemma 4.3.1, M 2 Gn M 2 Hn = o{Fn)- 
Since Hn = o{Fn), Lemma 4.3.2 shows that Ta^Hn = o{Fn)- Yet 
another application of Lemma 4.3.2 shows that 

Fu^Gn = 'y ^ F}j^Fn-k + l + o{Fn) 

l^k^n 

= Fn-k + l + o{Fn) ■ 
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This proves Theorem 2.3.2 when the weights are nonnegative. 

Removal of the sign restriction is accomplished via the decom¬ 
position of a weighted average into two weighted averages in which 
each sum has all weights with the same sign as was done at the end 
of the proof of Theorems 2.3.1 and 2.3.3. 


Acknowledgements. We thank Chenhua Zhang for shortening our 
original proof of Lemma 4.2.4. 
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